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Abstract 

CO adsorption on Cu(lll) and Cu(OOl) surfaces has been studied within ab-initio 
density functional theory (DFT). The structural, vibrational and thermodynamic 
properties of the adsorbate-substrate complex have been calculated. Calculations 
within the generalized gradient approximation (GGA) predict adsorption in the 
threefold hollow on Cu(lll) and in the bridge-site on Cu(OOl), instead of on-top as 
found experimentally. It is demonstrated that the correct site preference is achieved 
if the underestimation of the HOMO-LUMO gap of CO characteristic for DFT is 
correct by applying a molecular DFT-I-U approach. The DFT-I-U approach also 
produces good agreement with the experimentally measured adsorption energies, 
while introducing only small changes in the calculated geometrical and vibrational 
properties further improving agreement with experiment which is fair already at 
the GGA level. 
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1 Introduction 



CO adsorption on transition metals has been studied already for half a century 
and is in general considered as being rather well understood. [1]. The models 
for chemisorption based on density functional theory (DFT) are generally 
accepted [2,3] and it seems that there is nothing to investigate. However, 
there are some cases where DFT fails. The adsorption of a CO molecule on 
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the Cu(lll) surface is one of them and the question why theory predicts 
adsorption at low coverage to take place in the hollow sites while experiments 
find the top site to be preferred is hence of general interest. 

In a recent paper [4] we have tried to cover the state-of-the-art of the theoreti- 
cal description of CO adsorption on close-packed surfaces of transition and no- 
ble metals. Our calculated structural, vibrational and thermodynamical prop- 
erties have been compared to available experimental data. The dependence 
of these basic properties on the filling of the d-band as well as the correla- 
tion with the experimental site-preference have been discussed. Although the 
choice of the exchange-correlation functional leads only to a small difference in 
the calculated structural and vibrational properties, it influences very strongly 
the absolute value of the adsorption energy [4] . A large effort has been devoted 
in the past to study CO adsorption on Pt(lll) [5] where CO also prefers to 
sit in the on-top site at low coverage, but where standard ab-initio calcula- 
tions predict the hollow site again. Attempts to resolve this discrepancy have 
been undertaken by many researchers. Three ways to correct the prediction 
of the site preference have been proposed: (1) the molecular DFT-I-U method 
proposed by Kresse et al. [6] , (2) the use of hybrid functionals incorporating a 
fraction of exact (i.e. Hartree-Fock (HF)) exchange in the DFT functional [7], 
and (3) relativistic corrections [8]. The former two approaches are based on the 
fact that the interaction between the empty CO 2tt* orbitals and the metallic 
d-band is overestimated due to the tendency of DFT to underestimate the gap 
between the highest occupied molecular orbital (HOMO - 5a) and the lowest 
unoccupied molecular orbital (LUMO - 27r*) of the CO molecule. Therefore 
the position of the 2n* orbital is shifted to higher values by an additional 
on-site Coulomb interaction U added to the DFT Hamiltonian or a mixing 
of DFT- and HF-exchange also leading to an increased level splitting. In con- 
trast, the third approach concentrates on relativistic effects. Relativistic effects 
are small for the CO molecule and therefore changes in the (/-band of the Pt 
substrate are considered as essential for restoring agreement with experiment. 
It is therefore of interest to re-investigate the problem of CO adsorption on 
the much lighter metal Cu where the same discrepancy between theory and 
experiment has been found. 

The literature on CO adsorption on Cu(lll) and Cu(OOl) is not as extensive as 
for transition-metal surfaces like Pt or Pd. The adsorption of CO on Cu(lll) 
is well understood from the experimental point of view. At low coverage and 
T = 95 K the terminal (top) site is preferred, giving a (V3 X 73) LEED 
pattern. At higher coverages, bridge and hollow sites are also populated [9,10]. 
A contraction of the first layer appears at 0.44 ML, leading to a (1.5xl.5)R18° 
structure [11]. In contrast ab-initio theory proposes the hollow site as the 
favorite site, but the difference in the adsorption energies is quite small (< 0.1 
eV/molecule) [12,4]. On the other hand, the heat of adsorption hardly changes 
up to a coverage of 0.44 ML [13]. 
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On the Cu(OOl) surface a c(2x2) structure of the CO overlayer was observed 
of in low energy electron diffraction (LEED) studies [14,15,16,17]. The CO 
molecule is adsorbed in on-top, but the adsorption induces a negative shift of 
the workfunction at saturation coverage [18]. For CO adsorption on Cu(OOl) 
surfaces the theoretical cluster study by van Daelen et al. [19] reports the 
correct adsorption site, but the adsorption energies are too high (top site, local 
functional: Eads = —1-56 eV) or too low (top site nonlocal functional: Eads = 
—0.39 eV) compared to experiment. Lewis and Rappe [20] report for the 
same site (top) a lower adsorption energy, but still too high (in absolute value) 
compared to experiment (Eads — —1.239 eV). Although Ge and King [21,22] 
found that bonding at the bridge site is stronger by 80 meV than in the top site, 
and a barrier between the top and the bridge site which was not observed in 
a helium atom scattering experiment (HAS) [23,24]. Finally, we know of only 
one attempt to study CO on the Cu(OOl) surface including the generahzed 
gradient approximation (GGA) for the exchange-correlation functional [25]. 
The authors demonstrated that the revised Perdew-Burke-Erzernhof (RPBE) 
functional [26,27] (which has been adjusted to fit CO adsorption data) leads to 
better agreement with experiment, but the adsorption energy is still stronger 
by about 0.2 eV (or 35 %) than the experimental value. 

In this paper we present an investigation of the adsorption of CO on the 
copper (111) and (001) surfaces using the molecular DFT+U method pro- 
posed by Kresse et al. [6]. As copper is much lighter than platinum, rela- 
tivistic corrections are expected to be less important, while the correction 
of the HOMO-LUMO gap of CO should lead to the correct site-preference 
also in this problematic case. Section 2 reviews the basic DFT methodology 
and introduces the DFT+U method. Sections 3.1 and 3.2 summarize the de- 
scription of the structural, vibrational and thermodynamic properties of the 
adsorbate-substrate complex within the DFT. Section 4 apphes the molecular 
DFT-I-U method and demonstrates that it corrects the prediction of the site 
preference agreement and leads to adsorption energies in agreement with the 
experiments, while introducing only small changes in the geometric and vibra- 
tional properties. In the last section we discuss the advantages of the method 
and summarize all important points. 



2 Methodology 

2. 1 DFT calculations 

The calculations presented in this study are performed within a plane-wave 
based density functional framework. We have used the Vienna Ab-initio Simu- 
lation Package (VASP) [28,29] and employed projected augmented-wave (PAW) 
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potentials [30,31]. For the description of exchange and correlation the LSD 
functional proposed by Perdew and Zunger [32] is used, adding (non-local) 
generalized gradient corrections (GGA) of various flavor (PW91 according to 
Perdew et al. [26], RPBE the Perdew-Burke-Ernzernhof functional as modi- 
fied by Hammer et al. [26,27]). The equilibrium lattice constants of Cu are 
acu = 3.664 A (PW91), acu = 3.669 A (RPBE), compared to an experimental 
value of acu = 3.61 A. In spite of small difference between the PW91 and 
RPBE lattice constant we have used the lattice constant calculated by PW91 
for the calculations with RPBE. 

The substrates are modelled by four layers of metal separated by a vacuum 
layer of approximately the same thickness. The two uppermost substrate layers 
and the CO molecule are allowed to relax. A c(2x4) cell is used for the Cu(lll) 
surface and a p(2x2) cell for the Cu(OOl) surface, resulting in a coverage of 
0.25 monolayers. 

The self-consistent densities are calculated by iterative diagonalization of the 
Kohn-Sham Hamiltonian. The Fermi population of the Kohn-Sham states is 
kfiT = 0.2 eV and all the energies have been extrapolated to k^T ^0 eV. 
We have varied the k-point mesh from (4x3x1) to (8x6x1) for the (111) 
surface and from (4x4x1) to (8x8x1) for the (001) surface. The results here 
are presented for the most dense k-point mesh. For all further computational 
details we refer to our previous paper on CO adsorption on transition and 
noble metals [4] . 

2.2 The molecular DFT+U method 

In the molecular DFT+U method [6] the DFT functional is modified in such 
a way as to shift the LUMO to higher energies. The form of the total-energy 
functional is inspired by the conventional DFT-I-U method adding an on-site 
Coulomb repulsion to the DFT Hamiltonian in order to account for strong 
electronic correlation effects [33,34], 

Egga+u = Egga + ^ E E« - «). (1) 

^ (T=l j=l 

Here i?GGA and Egga+u are the standard DFT functional in the generalized 
gradient approximation and the modified density functional and cr is a spin 
index. The nf stands for the occupancies of the two LUMO orbitals for 

up and down spin. The occupancies are determined by defining two projec- 
tion operators which are 1 for the LUMO orbitals of the free CO molecule 
and for all other molecular eigenstates. Within the projector-augmented 
wave (PAW) approach implemented in VASP the projection operators defin- 
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ing these occupancies are identified with the PAW projections operators acting 
in the augmentation spheres around the carbon and oxygen atoms only. The 
occupancies n1 are hence determined taking the sum of 



(2) 



n,lm 



over all one-electron states "^"^ with fixed spin a. The Pn,im are the PAW 
projector functions for the states with quantum numbers Im on atom n as 
defined in Ref. [31]. The real coefficients a^i„j have to be chosen such that 
n"^ is one for a one-electron state belonging to the 27T* LUMO of a gas-phase 
CO molecule and zero for all other states. If the molecular axis is oriented 
along z, the 27r* molecular orbital consists of two orbitals 2n*, 27r*. Hence 
the coefficients Q;o,p^ and ac,p^ have to be chosen such that the corresponding 
occupancies are zero for the Itt and one for the 27r* molecular orbitals, i.e. 



where the ^s stand for the 27r* and ItTx one-electron states. These relations 

define two equations for the coefficients ao,p.c ^-^d ac,pa; , and similar relations 
hold for the Hy orbitals. For any further detail we refer to the paper by Kohler 
and Kresse [35]. 

It is important to emphasize that this ansatz leaves all DFT total energies 
unchanged if the occupancies are zero or unity. Hence the bond-length, the vi- 
brational cigcnfrequencies and even the excitation energies of the free molecule 
remain entirely unchanged. The important point is that despite the unmodified 
ground-state energy and structure, the adsorption energy decreases substan- 
tially, and that, as shown in detail below, the site preference is corrected. 

Kresse et al. [6] point out that with a modest value oi U — 0.75 eV, the 
DFT+U approach yields a single-particle 5(J — > 27r* gap in close agreement 
with the excitation energy for transferring a single electron from the 5a to 
the 2tt* state, while in the GGA without the Hubbard-like correction, the 
calculated excitation energy (not to be identified with the triplet to singlet 
excitation) is larger than the single-particle gap by 0.4 eV. This difference 
refiects the tendency of DFT to favor fractional orbital occupancies and hence 
to facilitate fractional electron transfer from the HOMO to the LUMO orbital. 
Other possibilities to correct this failure of DFT are hybrid functional incor- 
porating a fraction of exact exchange or self-interaction corrected functionals. 



(3) 
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Fig. 1. Schematic sketch of the molecular eigenstates of the CO molecule. The 
DFT+U technique shifts the LUMO orbitals to higher energies, but the energies 
of the occupied orbitals remain the same. The inset shows the variation of the 
HOMO-LUMO for free CO with increasing U. 

3 CO adsorption on Cu within DFT 

In this section we will briefly go trough our DFT results for CO adsorption 
on the Cu (111) and (001) surface. The structural, the vibrational and ther- 
modynamic properties are presented. 

3.1 Cu(lll) surface 

In Table 1 the structural, vibrational properties and the adsorption energies 
calculated using both the PW91 and RPBE functionals are presented. 

We discuss the PW91 results first. CO adsorption on the Cu(lll) surface 
follows the geometrical trends which depend on coordination [4]. The C-0 
bond is expanded compared to the gas-phase due to partial occupation of 27r* 
orbital, the expansion increases with the number of Cu atoms to which the 
molecule binds. The Cu-C distance varies in the interval from 1.86 to 2.05 
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Table 1. Calculated structural, vibrational and thermodynamical properties of CO adsorbed in the high symmetry sites on Cu(lll) 
for a coverage of 9 = | ML (c(2x4)): Ea^s - adsorption energy, dc-o - carbon-oxygen bond length , hco - height of CO molecule 

above the surface, dcu-C - Cu-C bond length, A(ii2, A(i23 - change of the average inter-laycr spacing, 61,62 - buckling of 1** and 2"''^ 
layer, uc-O - intramolecular stretching frequency, vs-CO - surfacc-adsorbatc stretching frequency, A<I> - adsorption-induced change of 
the work-function. The results have been obtained with the PW91 functional, selected results calculated with the RPBE functional and 
PW91-I-U method (U = 1.25 eV) are given in rectangular and round parentheses. Experimental data taken from Refs. [36,37,9,38,10,39,40] 
are listed where available. 



A and increases with coordination. The buckhng of the first layer and the 
inward-relaxation of the top Cu-layer tend to decrease with coordination. 

The experimental stretching frequency of the CO molecule {pc-o) changes 
upon adsorption from the gas-phase value of vc-o — 2145 cm~^ to about 2070 
cm~^ ( red-shift Au ~ -75 cm~^), depending on coverage and temperature. 
The calculated stretching frequency of gas-phase CO frequency (2136 cm~^, 
Ref. [4]) shifts to 2046 cm~^ [Ai> ~ -90 cm~^). In spite of the difference of 15 
cm~^ between the measured and calculated red-shifts, the agreement between 
experiment and theory is very satisfactory, because the stretching frequency 
of adsorbed CO decreases as the temperature decreases and therefore the 
shift measured in room-temperature experiments is smaller. The same trend 
(shift downwards) is observed for increasing coverage [10,13]. The difference 
between the experimental and theoretical vm-co is also quite small ( ~ 10 
cm"^). Bridge and hollow adsorption lead to far larger redshifts of both the 
C-0 and the substrate-CO stretching frequencies. 

The workfunction of the clean surface is calculated to be $dcan^^^= 4.78 cV 
for PW91 which is lower than experimentally estimated 4.98 eV [40]. On-top 
adsorption of CO leaves the work function almost unchanged (A$ = 0.06 eV), 
while the stronger interaction in the sites with higher coordination leads to a 
pronounced increase of $ by 0.6 eV (bridge) and 0.7 eV (hollow). 

Altogether we note good agreement of the vibrational frequencies calculated 
for a linear adsorption geometry with experiment, while the too large adsorp- 
tion energy and the too short Cu-C distance demonstrate that the ovcrbinding 
tendency of the LDA is not completely removed by the PW91 form of the GGA 
functional. The most serious discrepancy, however, is the preference for bind- 
ing in the hollow sites, the difference in the adsorption energy compared to the 
top site is 0.10 eV/molecule. Hammer et al. [27] have noted that a shght modi- 
fication of the PBE functional [26] (which does not affect any of the sum rules) 
largely improves the prediction of the adsorption energies. This is confirmed 
by our calculations. Using the RPBE functionals we find adsorption energies in 
very good agreement with experiment (see Table 1). The energetic preference 
for adsorption in the fee hollow is reduced to 0.03 eV/molecule, but still the 
wrong adsorption site is predicted. The change in the GGA functional has only 
a minor influence on the geometric and vibrational properties: the elongation 
of the C-0 bond length is further increased by 0.006 A, the distance between 
the C-atom and the nearest Cu atom by 0.01 A (top) to 0.03 A (hollow). The 
most important change occurs for the outermost interlayer distance of the Cu 
substrate which is expanded by about 3 % compared to the PW91 results. The 
C-0 stretching frequencies are reduced on average by about 10 cm~^ for the 
gas-phase and adsorbed CO, so that the red-shift remains the same. A slightly 
more pronounced softening is observed for the adsorbate-substrate stretching 
frequency, but altogether the change of the GGA functional has only a small 
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influence on geometry and dynamics. 



3.2 Cu(OOl) surface 

For CO adsorption on tiiis surface tlie number of tlieoretical investigations 
exceeds tfiat of experimental studies - which could indicate that there are still 
unresolved questions. Our study covers adsorption in the three high-symmetry 
surface sites and the structural, vibrational and thermodynamic properties of 
CO on Cu(OOl) are compiled in Table 2. 

If we compare CO adsorption on Cu(lll) and Cu(OOl) similar trends are ob- 
served: the elongation of the C-0 bond and the Cu-C bond-lengths depend 
on coordination. With increasing coordination the adsorbatc-substrate bond 
length increases, the inward relaxation of the top layer and the buckling am- 
phtude are reduced. Again we find a pronounced decrease of the adsorption 
energies if the PW91 functional is replaced by the RPBE functional. The 
reduction is stronger for the bridge and hollow sites, so that the RPBE cal- 
culation shows even a very slight preference for the experimentally observed 
on-top geometry by 0.02 eV compared to the bridge site. This energy differ- 
ence is, however, much smaller that the 80 meV determined experimentally 
by Ge and King [21,22]. 

The adsorption geometries obtained with the different exchange correlation 
functionals differ only very slightly. The C-0 bond length in the on-top con- 
figuration is 1.156 A for the PW91 functional, and 1.162 A for the RPBE 
functional. Similarly the metal-carbon distance is expanded from 1.85 A with 
PW91 to 1.87 A with RPBE. Agreement is slightly improved with the RPBE 
functional, but the difference is smaller than the experimental uncertainty. 

The calculated red-shifts of the CO stretching mode are slightly smaller on 
the Cu(OOl) than on the Cu(lll) surface for on-top {Au ~ —88 cm~^) and 
bridge (Az/ ~ —210 cm^^) adsorption. Again the calculation overestimates 
the redshift (experiment: Au ~ —66 cm~^ for on-top CO). For adsorption in a 
hollow, however, one finds an increased red-shift of Av ~ -418 cm~^, compared 
to Av ~ -286 cm~^ on Cu(lll), reflecting the smaller adsorbate height in a 
four-fold than in a thrcc-fold hollow and the larger stretch of the CO bond. 
The small changes of the frequencies calculated with the RPBE instead of the 
PW91 functional are of the same order as on the Cu(lll) surface. 

With $2n°^^ = 4.57 eV (PW91) and 4.49 eV (RPBE) the workfunction of the 
clean Cu(OOl) surface is smaller than for the Cu(lll) surface. Experimentally 
values covering the interval from 4.59 to 4.83 eV have been reported [18,42,43]. 
Upon adsorption the workfunction increases quite strongly, depending on co- 
ordination of the adsorption site: using PW91 we find A$ = 0.29 eV (top). 
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Table 2. Calculated structural, vibrational and thermodynamical properties of CO adsorbed in the high symmetry sites on Cu(OOl) for 
coverage © = | ML. The nomenclature is the same as in Table 1. Results obtained using the both the PW91, the RPBE (in rectangular 
parentheses) and the PW91+U (in round parentheses, U = 1.25 eV) functionals are listed. The experimental data are taken from Refs. 
[36,16,14,15,41]. 



0.65 eV (bridge), and 0.89 eV (hollow), RPBE shifts are slightly larger by 
about 0.05 eV. 



4 Improving site preference within DFT+U method 




Fig. 2. The adsorption energies of CO on the Cu(lll) surface as a function of the 
parameter U for the molecular DFT+U method. 

The DFT+U method enables us to study CO adsorption with different amounts 
of interaction between the unoccupied 27r* molecular orbital and the metal 
d-band. We have studied the variation of the adsorption energy with an in- 
creasing value of the on-site potential U. Fig. 2 demonstrates the change in 
the adsorption energy of CO on the Cu(lll) surface with increasing U. The 
calculations are based on the PW91 functionals - if the DFT+U approach 
is based on the RPBE functional, the adsorption energies are drastically un- 
derestimated. The adsorption energy linearly decreases with increasing U and 
the rate of reduction is greater for higher coordinated sites. Site-preference 
between top and bridge is reversed already around U=0.05 eV and between 
top and hollow around U=0.47 eV. For U between ~ 0.9 eV and ~ 1.25 eV 
the calculated adsorption energies fall into the interval covered by the experi- 
mental values. The same analysis for CO/Cu(001) is given in Fig. 3. Stability 
of on-top adsorption over the bridge site is achieved at U> 0.15 eV, to match 
the experimentally determined difference in the adsorption energies (^/SEads ~ 
80 meV), a value of U~0.85 eV is required, agreement with the experimen- 
tally measured adsorption energies is reached for U ~ 1.25 eV. It is certainly 
encouraging the the same values of the on-site potential fit the CO adsorption 
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data on both surfaces and that the optimal value of U determined from the 
adsorption study agrees reasonably well with that needed to bring the one- 
electron HOMO-LUMO gap in agreement with the excitation energies (cf. 
Sec. 2). A further confirmation that the molecular DFT+U approach grasps 
the correct physical mechanism comes from the observation that the difference 
of the on-top adsorption energies on the Cu(lll) and Cu(OOl) surfaces is in- 
dependent of U (see Fig. 2). For adsorption in the fourfold hollow on Cu(OOl), 
the decrease of the adsorption energy with U is faster than for adsorption in 
one of the threefold hollows on Cu(lll), demonstrating the dependence of the 
27r* backdonation on the coordination of the adsorbate. 



The structural and vibrational properties of CO adsorbed in the top site on 
Cu(lll) and Cu(OOl) calculated by the DFT+U approach are listed in Ta- 
bles 1 and 2. For these calculations we have used a Hubbard U parameter 
equal to 1.25 eV which is the highest acceptable value. The C-0 bond (1.156 
A calculated with PW91) is slightly shorter (1.152 A) and the Cu-C bond 
(1.85 A) is elongated to (1.87 A) for CO adsorption on Cu(lll). The Cu-C 
bond length calculated by DFT+U has the same value as the calculation with 
RPBE functional (1.87 A). The structure of CO on Cu(OOl) is influenced in 
the same manner. As the DFT+U corrections affect only the adsorbate, it 
is not surprising that the calculated interlayer relaxations are the same as 
predicted using PW91 alone. 
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Fig. 3. The adsorption energies of CO on the Cu(OOl) surface as a function of the 
parameter U in the molecular DFT+U method. 

The stretching frequency of CO molecule {uc-o) calculated using the DFT+U 
approach is shifted to higher frequencies by 20-40 cm~^. The DFT+U frequen- 
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cies of uq-o arc closer to the experimental values. The stretching frequency 
of CO (lyc-o) ii^i top on Cu(lll) is 2075 cm^^ and the corresponding exper- 
imental values are 2071 and 2077 cm~^. Similarly, the calculated /Uc-o in 
top on Cu(OOl) is 2078 cm~^ and the experimentally observed values are 
2079 and 2085 cm~^. The stretching frequency of the adsorbate-substrate 
bond, (us^co), is decreased by ~ 10-20 cm~^ compared to PW91. Hence, 
the DFT+U approach leads to comparable results as the calculations with 
RPBE, from 20 to 35 cm~^ below the experimental values. 

A stronger influence of the corrected HOMO-LUMO gap is found for the 
workfunction. The shorter C-0 bond resulting from the DFT+U calculations 
decreases the workfunction compared to a PW91 calculation due to the reduc- 
tion of the dipole moment by ~ 0.2-0.3 eV. The workfunction change of the 
system where the CO molecule sits on top of Cu becomes negative ( A$ =-0.22 
eV) on the Cu(lll) surface and shghtly positive (A$ ~ 30 meV) for on-top 
CO on Cu(OOl). 

We can conclude that the DFT-I-U method slightly modiflcates the structural 
and vibrational properties which are in better agreement with experiment than 
the values obtained by PW91 and RPBE exchange-correlation functional. 



5 Conclusion 

A detailed analysis of the DFT predictions for CO adsorption on the Cu(lll) 
and Cu(OOl) surface was presented. Calculations with two different GGA func- 
tional produce rather accurate adsorption geometries and vibrational eigen- 
frequencies for CO adsorbed in the on-top position. However, using the PW91 
functional too large adsorption energies and a wrong site preference are pre- 
dicted on both surfaces. Replacing the PW91 by the RPBE functionals hardly 
changes the vibrational frequencies and the geometry (except the relaxation 
of the top layer) and results in adsorption energies in much better agree- 
ment with experiment. However, the site-preference is corrected only for the 
Cu(OOl), but not for the Cu(lll) surface. 

It was demonstrated that the site preference may be corrected using a molec- 
ular DFT+U approach which increases HOMO-LUMO separation in the CO 
molecule. This leads to a reduced electron donation into the antibonding 27r* 
orbital and reduces the strength of the chemisorption, the reduction being 
more pronounced for a higher coordination. The weaker adsorption strength 
results also in a higher stretching frequency of the adsorbed molecule, also 
improving agreement with experiment. Altogether we find that the molecular 
DFT+U method is a very convenient method to correct for the underestima- 
tion of the HOMO-LUMO gap characteristic for the DFT and for achieving a 
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better description of the adsorption of difficult molecules such as CO. 
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